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1.  Introduction 


Improving  delivery  accuracy  of  squad-level  weapons  is  the  motivation  for  this  report.  Guiding 
projectiles  is  a  means  of  compensating  for  delivery  errors  due  to  effects  such  as  wind,  weapon 
aiming,  or  muzzle  velocity  variation.  Maneuvering  projectile  flight  is  a  critical  element  of  any 
guidance  technology.  The  focus  of  this  report  is  maneuvering  small-diameter  projectiles. 

Guns  are  often  rifled  to  impart  spin  to  projectiles,  which  results  in  gyroscopic  stability  (7,  2). 

One  advantage  of  spin-stabilized  projectiles  over  fin-stabilized  projectiles  is  reduced  drag.  Flight 
stability  of  controlled,  spin-stabilized  projectiles,  however,  is  more  complex  than  ballistic  flight 
stability  (3,  4).  Maneuvers  directed  perpendicular  to  the  gravity  vector  increase  the  chance  for 
unstable  flight.  Additionally,  high-spin  rate  complicates  the  actuation  technology  (e.g.,  higher 
bandwidth  requirements).  Portions  of  the  projectile  or  the  control  mechanism  may  be  decoupled 
from  the  main  projectile  body  to  mitigate  these  issues  (5-9).  In  contrast  to  fin-stabilized 
projectiles,  the  optimal  placement  for  maneuverability  of  spin-stabilized  projectiles  is  at  the  rear 
of  the  body  due  to  the  center-of-pressure  location  (10, 11). 

The  launch  velocity  of  man  portable  systems  is  often  limited  by  weapon  recoil.  There  is  a 
maximum  acceptable  impulse  of  the  launch  event  on  the  shoulder  of  the  Soldier.  Low-launch 
velocity  translates  to  low-dynamic  pressure  available  for  maneuvering  an  aerodynamic  control 
mechanism.  Thrusters,  nonaerodynamic  in  nature,  are  viable  options  in  low-dynamic  pressure 
flight  regimes  (12-1 7).  The  angular  motion  of  the  projectile  induced  by  impulsive  maneuvers 
may  be  deleterious  to  seeker-based  guidance  performance.  Current  flow  control  technologies 
feature  limited  control  authority  (18). 

The  size,  weight,  power,  and  cost  requirements  for  the  maneuver  technologies  in  small-diameter 
projectiles  are  challenging.  Man  portable  weapons,  from  assault  rifles  to  grenade  launchers, 
shoulder  fired  weapons,  and  mortars  are  often  about  80  mm  in  diameter  or  smaller.  Actuators, 
sensors,  electronics,  embedded  processors,  and  a  power  source  must  fit  along  with  other  potential 
subsystems  such  as  fuzing  and  warhead.  Current  fin  or  canard  actuation  technologies  are  difficult 
to  package  in  less  than  80-mm-diameter  projectiles  (19-22). 

The  high-structural  loads  applied  to  these  components  during  gun  launch  impose  severe  constraints 
(23, 24).  Reducing  the  size  and  number  of  moving  parts,  as  well  as  properly  supporting  components, 
is  often  necessary  for  ensuring  survivability  and  reliability  in  the  gun-launched  environment. 

The  goal  of  this  report  is  to  explore  the  aeromechanics  of  a  small  guided  projectile  concept  to 
assess  the  efficacy  in  improving  delivery  accuracy.  A  rotating  wing  actuator  is  the  basis  of  this 
new  maneuver  concept.  The  novel  contributions  are  the  flight  concept,  theoretical  modeling  of 
the  concept,  computational  and  experimental  aerodynamic  characterization,  and  assessment  of 
control  authority  with  respect  to  ballistic  dispersion. 
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This  report  first  outlines  the  maneuver  concept  and  geometry  of  the  projectile  and  control 
parameters.  Next,  aerodynamic  data  obtained  through  experimental  and  computational 
techniques  are  provided.  Aerodynamics  were  used  in  a  model  as  input  to  the  flight  mechanics. 
Lastly,  launch  and  flight  system  simulations  are  shown  which  quantify  the  ballistic  dispersion, 
guided  flight  stability  and  control  authority. 


2.  Concept 


The  maneuver  concept,  shown  in  the  figure  1,  applies  to  spin-stabilized  projectiles.  A  wing¬ 
shaped  actuator  aligned  with  the  spin  axis  of  the  projectile  is  used  to  constantly  rotate  an 
aerodynamic  surface,  such  as  a  miniature  wing,  with  respect  to  the  projectile  body.  Over  one 
revolution  of  projectile  spin  the  actuator  rotates  the  aerodynamic  surface  or  wing  from  exposed 
to  the  airstream  (effecting  lateral  maneuver  as  shown  in  figure  1)  to  stowed  internally  within  the 
projectile  body.  Rotating  the  aerodynamic  surface  opposite  the  projectile  spin  achieves  a 
consistent  maneuver  direction.  Phase  shifting  the  wing  rotation  with  respect  to  the  projectile 
rotation  changes  the  maneuver  direction. 


Figure  1.  Rotating  wing  maneuver  concept. 

Further  pictorial  description  of  the  maneuver  concept  is  provided  in  figure  2.  A  complete 
revolution  of  spin  advances  in  45-deg  roll-angle  increments  in  a  counter-clockwise  manner 
around  the  figure.  A  two-actuator  variant,  viewed  from  the  base  of  the  projectile  in  an  Earth- 
fixed  frame,  is  illustrated.  The  projectile  rotates  in  the  clockwise  direction  and  the  actuation 
assembly  rotates  in  the  counter-clockwise  direction.  Each  successive  rendering  illustrates  the 
manner  in  which  the  rotational  actuator  exposes  the  wing  to  the  airstream  to  produce  a  lateral 
maneuver. 
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Figure  2.  Snapshots  (viewed  from  base  in  earth-fixed  reference  frame)  over  one  revolution  of  rotating  wing 
maneuver  concept  with  two  actuators. 

This  maneuver  concept  contrasts  with  conventional  aerodynamic  surface  actuation  technologies 
wherein  rotational  actuators  travel  less  than  one  full  rotation  (e.g.,  servomechanisms)  to  deflect 
aerodynamic  surfaces.  Linear  actuation  techniques  are  not  as  attractive  in  the  high-spin  rate 
environment  due  the  power  requirements  (25,  26).  The  approach  outlined  in  this  report  marries 
low-cost  actuation  technology  with  spin-stabilized  flight  dynamics  for  a  novel  solution. 
Rotational  actuators  with  spin  rates  in  the  10-1000-Hz  range  are  commercially  available  in  high 
volume  at  low  cost  with  small  size,  weight,  and  power  characteristics.  Companion  efforts  are 
investigating  the  mechatronics,  control,  and  performance  of  the  maneuver  technology  (27). 
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3.  Aerodynamic  Characterization 


The  projectile  geometry  is  presented  in  figure  3.  A  hemispheric  nose,  adequate  for  subsonic 
flight,  is  followed  by  a  2.25-caliber  (cal.)  cylindrical  section  prior  to  a  7-deg  boattail  for  wake 
drag  reduction.  The  engraving  band  was  0.125  cal.  long  and  0.02875  cal.  in  height.  The  specifics 
of  this  configuration  were  chosen  based  on  an  analysis  of  subsystem  requirements  to  maximize 
range  considering  the  following  launch  and  flight  characteristics:  caliber,  projectile  length- to- 
diameter  ratio,  stabilization  method,  recoil  impulse,  and  propulsion  (gun  or  gun  and  rocket)  (28). 
A  center-of-gravity  location  1.675  cal.  from  the  nose  was  used  for  all  calculations. 


A  parametric  analysis  of  the  wing  was  conducted  to  define  optimal  characteristics.  The 
parameters  under  investigation  were  the  chord  (C),  span  (S),  deflection  angle  (5),  airfoil  profile, 
and  stand-off  from  the  projectile  body  (H).  The  geometry  for  the  wing  is  provided  in  figure  4. 
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A  comprehensive  aerodynamic  characterization  of  the  projectile  and  wing  were  undertaken.  The 
approach  was  to  initially  perform  semiempirical  aeroprediction  (SEAP)  in  the  Projectile, 
Rockets,  and  Ordnance  Design  and  Analysis  System  (PROD AS)  Spinner  code  on  the  projectile 
body  (29).  Projectile  and  wing  static  aerodynamic  data  were  collected  during  wind  tunnel  (WT) 
experiments.  These  data  served  to  identify  optimal  wing  characteristics  and  validate 
computational  fluid  dynamics  (CFD)  predictions.  Modeling  down-selected  geometries  and 
performing  CFD  provided  dynamic  derivatives  and,  along  with  the  WT  data,  enabled  the 
aerodynamic  model  to  be  formulated.  Characterizations  were  performed  of  the  projectile  body 
and  of  the  projectile  body  with  the  wing.  To  isolate  the  maneuver  aerodynamics,  projectile  body 
aerodynamic  data  were  subtracted  from  the  projectile  body  with  wing  data.  Thus,  wing  and  any 
wing-body  interference  effects  were  modeled  as  part  of  the  maneuver  aerodynamics. 


The  body-fixed  coordinate  system,  defined  in  the  figure  5,  was  used  to  quantify  aerodynamics. 
The  coordinate  axes,  and  subsequently  the  velocity,  forces,  and  moments,  follows  the  right-hand 
rule.  Aerodynamic  angles  (pitch  angle-of-attack,  a,  yaw  angle-of-attack,  /?,  total  angle-of-attack, 
a,  and  aerodynamic  roll  angle,  <pA)  are  computed  based  on  the  body-fixed  velocity  components 
([u  v  w]T). 


P 

a 


a  =  asm 


w 


Vu2  + 


w' 


=  asm 


=  asm 


v 


-^lu2  +  V2  +  w2 


Vu2  + 


w' 


Vlt2  +  V2  +  w2 


4>a 


=  atari 


(1) 
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Subsonic  experiments  were  conducted  in  a  blow-down  WT  with  a  28  in  high  x  40  in  wide  test 
section  housed  at  the  Edgewood  Area  of  Aberdeen  Proving  Ground,  MD.  Models  were 
constructed  with  rapid-prototyping  tools.  The  projectile  diameter  was  70  mm  and  was  made  from 
polycarbonate  with  fused  deposition  modeling.  The  wings  were  made  from  Somos  Nanotool 
with  stereolithography.  The  models  were  mounted  on  a  sting  with  a  3/8-in-diameter  five 
component  balance  (Modem  Machine  &  Tool  Co.).  The  balance  featured  sensitivities  of  10  lb 
for  axial  force,  8  lb  for  normal  force,  5  lb  for  side  force,  and  8  in-lb  for  pitching  moment  and  side 
moment.  An  image  of  the  model  with  balance  on  the  sting  in  the  tunnel  is  shown  in  figure  6. 
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Figure  6.  Instrumented  model  mounted  in  WT  facility  test 
section. 

Data  were  obtained  in  the  WT  at  angles-of-attack  of  ±12  deg  in  1-  or  2-deg  increments.  The 
Mach  number  was  approximately  0.2  since  this  is  of  tactical  interest  for  low-velocity  grenade 
weapons  (30).  The  WT  key  and  matrix  is  given  in  tables  1  and  2.  Two  bodies  were  examined;  the 
first  was  a  replica  of  the  projectile  geometry  outlined  in  figure  3  and  the  second  was  a  similar 
body  with  a  cavity  to  accommodate  the  wing.  The  effect  of  airfoil  camber  was  investigated  along 
with  planform.  The  wing  deflection  angle  was  studied  to  assess  stall.  Standoff  of  the  wing  from 
the  projectile  body  was  varied.  The  effectiveness  of  the  wing  as  a  function  of  the  amount  of  wing 
exposed  to  the  airstream,  or  stowed  roll  angle  (0,  measure  of  amount  of  control  surface  exposed 
to  airstream)  was  obtained.  The  stowed  roll  angle  can  be  visualized  in  figure  7.  Here,  the 
decrease  in  the  wing  area  exposed  to  the  airstream  as  stowed  roll  angle  increases  from  0  to 
90  deg  in  30-deg  increments  can  be  seen  with  the  model  in  the  WT.  When  describing  a  particular 
configuration  the  format  will  be  to  concatenate  the  identification  number  for  successive 
parameters  in  the  WT  key.  For  example,  configuration  1 1 120  is  the  projectile  with  cavity, 
NACA*  4415  airfoil,  0.48-cal.  chord  x  0.61 -cal.  span,  and  8-deg  deflection.  Data  were  also 
collected  at  aerodynamic  roll  angles  of  0,  22.5,  45,  67.5,  and  90  deg. 


* 

National  Advisory  Committee  for  Aeronautics. 
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Table  1.  WT  matrix  key. 


Parameter 

Identification 

No. 

Description 

Body 

0 

Projectile  body 

1 

Projectile  body  with  cavity 

Airfoil 

0 

NACA  0015 

1 

NACA  4415 

Planform 

0 

0.61 -cal.  chord  x  0.61 -cal.  span 

1 

0.48-cal.  chord  x  0.61 -cal.  span 

2 

0.32-cal.  chord  x  0.61 -cal.  span 

3 

0.48-cal.  chord  x  0.48-cal.  span 

4 

0.32-cal.  chord  x  0.48-cal.  span 

Deflection 

0 

4  deg 

1 

6  deg 

2 

8  deg 

3 

10  deg 

Standoff 

0 

5.28  mm  from  body  outer  diameter 

1 

0.3  mm  from  body  outer  diameter 

Stowed  roll  angle 

0 

0  deg 

1 

30  deg 

2 

60  deg 

3 

90  deg 

Aerodynamic 
roll  angle 

0 

0  deg 

1 

22.5  deg 

2 

45  deg 

3 

67.5  deg 

4 

90  deg 

5 

60  deg 

6 

30  deg 

Table  2.  WT  matrix. 


Body 

Airfoil 

Planform 

Deflection 

Standoff 

Stowed  Roll 

Aero  Roll 

0 

NA 

NA 

NA 

NA 

NA 

NA 

1 

NA 

NA 

NA 

NA 

NA 

0 

1 

1 

1 

2 

0 

0 

0 

1 

0 

1 

2 

0 

0 

0 

1 

1 

2 

2 

0 

0 

0 

1 

1 

3 

2 

0 

0 

0 

1 

1 

4 

2 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 

3 

0 

0 

0 

1 

1 

1 

2 

1 

0 

0 

1 

1 

1 

2 

0 

0 

1 

1 

1 

1 

2 

0 

0 

2 

1 

1 

1 

2 

0 

0 

3 

1 

1 

1 

2 

0 

0 

4 

1 

NA 

NA 

NA 

NA 

NA 

4 

1 

1 

1 

2 

0 

1 

5 

1 

1 

1 

2 

0 

2 

6 

1 

1 

1 

2 

0 

3 

0 

0 

NA 

NA 

NA 

NA 
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Figure  7.  Stowed  roll  angle  (left-to-right:  0,  30,  60,  and  90  deg). 


Data  collected  at  discrete  angles-of-attack  were  inspected  and  least-squares  fit  with  an 
appropriate  polynomial  expansion  in  angle-of-attack.  Some  examples  of  WT  data  and  subsequent 
polynomial  fits  are  shown.  The  fits  for  body-wing  and  body  alone  were  subtracted  as  described 
earlier  to  isolate  the  control  effects.  Figure  8  shows  the  normal  force  coefficient  for  different 
airfoils.  The  projectile  body  data  and  fit  are  in  black.  The  NACA  0012  data  are  in  blue  and  the 
NACA  4412  data  are  purple.  Points  represent  WT  data  and  curves  represent  fits.  Solid  lines 
represent  the  control  force  and  moment  obtained  by  subtracting  the  body-wing  from  the  body- 
alone  data.  As  expected,  the  body-alone  data  goes  through  the  origin.  When  adding  the 
aerodynamic  asymmetry  of  the  deflected  wing,  however,  a  trim  force  and  change  in  nonlinearity 
with  angle-of-attack  emerges.  These  data  show  that  adding  camber  increases  the  effectiveness  of 
the  control  by  almost  a  factor  of  2  in  the  trim  force. 
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The  pitching  moment  for  various  airfoils  is  given  in  figure  9.  Again,  the  body-alone  data  goes 
through  the  origin.  The  slope  of  the  body  data  is  positive  since  the  projectile  is  statically 
unstable.  Addition  of  the  wing  greatly  changes  the  nonlinearity  and  again  features  a  nonzero 
intercept.  Pitching  moment  for  the  NACA  4412  indicates  a  greater  trim  moment,  which  results  in 
more  control  authority. 


Figure  9.  Pitching  moment  coefficient  for  different  airfoils. 


Figures  10  and  1 1  provide  the  normal  force  and  pitching  moment  for  the  various  deflection 
angles.  Trends  are  similar  to  those  described  for  the  airfoil  parameter.  Addition  of  the  wing 
induces  a  trim  force  and  moment  and  change  to  the  nonlinearity  with  angle-of-attack.  The  wing 
effectiveness  increases  with  deflection  angle  since  the  lift  on  a  wing,  scales  with  effective  angle- 
of-attack  (projectile  body  and  deflection  angle  combined).  Stall  on  the  wing  does  not  appear 
significant  even  for  effective  angles-of- attack  greater  than  12  deg. 
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Figure  10.  Normal  force  coefficient  for  different  deflections. 
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Figure  11.  Pitching  moment  coefficient  for  different  deflections. 
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The  aerodynamic  data  for  all  configurations  from  WT  experiments  is  compiled  in  table  3.  In 
general,  the  wing  increases  drag  while  also  providing  the  necessary  normal  force  for  maneuver. 
The  cambered  airfoil  outperforms  the  airfoil  without  camber.  Larger  planform  increases 
effectiveness.  Control  effectiveness  increased  with  deflection,  although  the  results  at  10  deg  were 
not  appreciably  higher  than  those  at  8  deg.  The  larger  standoff  distance  was  critical  to  getting  the 
flow  around  the  wing  sufficiently  separate  from  the  flow  near  the  projectile. 

Table  3.  Aerodynamics  from  WT  for  all  configurations. 


Configuration 

Cx0 

Cx  2 

Cvo 

cNa 

CN  3 
or3 

r 

r 

u ma 

r 

ma 3 

0 

0.169782 

1.011851 

0.010264 

1.829601 

-4.94869 

0.005764 

1.932517 

-1.99989 

11120 

0.070048 

0.2794 

0.172031 

0.437611 

-0.75652 

-0.08572 

-0.10841 

-0.44945 

10120 

0.03475 

0.474122 

0.106159 

0.35477 

2.114011 

-0.03716 

-0.11304 

0.665462 

11220 

0.030063 

0.393898 

0.112284 

0.311331 

-0.09466 

-0.05199 

-0.06694 

-0.21999 

11320 

0.036989 

0.367927 

0.120833 

0.199218 

3.137703 

-0.05406 

-0.10749 

2.116242 

11420 

0.019322 

0.26236 

0.071913 

0.011366 

4.718838 

-0.02909 

-0.07454 

2.58084 

11100 

0.042328 

0.386538 

0.129257 

0.304703 

3.014388 

-0.05265 

-0.17627 

2.30084 

11110 

0.054294 

0.155404 

0.14464 

0.36058 

0.283997 

-0.06663 

-0.17575 

1.483676 

11130 

0.075184 

0.223102 

0.186378 

0.415412 

0.537722 

-0.09615 

-0.18314 

0.603045 

11121 

0.042328 

0.386538 

0.129257 

0.304703 

3.014388 

-0.05265 

-0.17627 

2.30084 

The  effect  of  stowed  roll  angle  is  demonstrated  in  table  4.  As  less  of  the  wing  is  exposed  to  the 
airstream,  the  drag,  lift,  and  moment  decrease.  The  wing  effectiveness  does  not  fall  sharply  with 
stowed  roll  angle;  however,  since  this  concept  is  able  to  keep  appreciable  wing  in  the  freestream 
even  at  higher  stowed  roll  angles  as  seen  in  figure  7.  A  roll  torque  is  expected  to  be  imparted  to 
the  projectile  during  the  partially  stowed  stages  of  the  cycle;  however,  roll  moments  were  not 
measured  in  the  WT.  Roll  moments  were  evaluated  via  CFD  and  discussed  in  the  following. 

Table  4.  Aerodynamics  from  WT  for  configuration  1 1 120  for  all  stowed  roll  angles. 


Cx0 

cv 

CV0 

CNa 

CN  3 
or3 

r 

Lmo 

r 

r 

ma 3 

0  deg 

0.070048 

0.2794 

0.172031 

0.437611 

-0.75652 

-0.08572 

-0.10841 

-0.44945 

30  deg 

0.059063 

0.038751 

0.151855 

0.427696 

0.636527 

-0.09272 

-0.22951 

0.171845 

60  deg 

0.026487 

-0.0745 

0.07989 

0.436496 

2.563468 

-0.06165 

-0.20457 

0.082161 

90  deg 

0.024284 

-0.19 

0.055608 

0.299621 

5.386849 

-0.04791 

-0.30041 

0.262917 

The  CFD  aerodynamics  characterization  (31)  consisted  of  a  WT  validation  series  and  a  down- 
selected  concept  characterization  series.  Configuration  11120  was  the  down-selected  concept  for 
the  remainder  of  the  CFD  analysis  due  to  its  attractive  characteristics  based  on  the  WT 
experiments. 

The  CFD  characterization  was  performed  using  the  commercially  available  CFD++  code,  version 
12.1  (32).  CFD++  solves  the  three-dimensional,  compressible  Reynolds- Averaged  Navier-Stokes 
(RANS)  equations  using  a  finite  volume  method.  The  inviscid  flux  function  was  a  second-order, 
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upwind  scheme  using  a  Harten-Lax-Van  Leer-Contact  (HLLC)  Riemann  solver  and  a 
multidimensional  Total-Variation-Diminishing  (TVD)  continuous  flux  limiter  (32).  The  solution 
is  advanced  to  steady-state  using  a  point-implicit  time  integration  scheme  with  local  time 
stepping  and  multigrid  convergence  acceleration.  The  low-speed  preconditioning  version  of  the 
compressible  solver  was  used,  as  it  is  appropriate  for  Mach  numbers  less  than  about  0.3.  CFD++ 
has  many  common  turbulence  models  available.  Menter’s  Shear  Stress  Transport  (SST)  model 
(33)  was  chosen  for  this  study  based  on  some  initial  comparisons  with  predictions  using  the 
realizable  k-s  model  (34)  and  the  WT  data  for  one  body-alone  configuration.  Simulations  were 
performed  primarily  on  the  SGI  ICE  X  System  (SPIRIT)  located  at  the  U.S.  Air  Force  Research 
Laboratory  (AFRL)  Department  of  Defense  (DOD)  Supercomputing  Resource  Center  (DSRC)  at 
Wright-Patterson  Air  Force  Base,  OH.  Simulations  were  also  undertaken  on  the  Cray  XE6 
(GARNET)  located  at  the  U.S.  Army  Engineering  Research  and  Development  Center  (ERDC) 
DSRC  in  Vicksburg,  MS,  and  the  IBM  iDataPlex  (PERSHING)  located  at  the  U.S.  Army 
Research  Laboratory  (ARL)  DSRC  at  Aberdeen  Proving  Ground,  MD. 

The  first  phase  of  the  CFD  characterization  was  a  validation  series  consisting  of  selected 
configurations  tested  in  the  WT.  The  Mach  number  was  0.2  and  the  WT  model  was  1.75  scale 
(70  mm).  Figure  12  shows  several  of  the  configurations  that  were  simulated.  CFD 
characterization  was  performed  for  the  solid  body  alone,  the  body  alone  with  empty  actuator 
cavity,  and  the  body  with  actuator  wing  deflections  of  4,  6,  8,  and  10  deg.  CFD  predictions  for 
configuration  1 1 120  actuator  wing  were  performed  for  aerodynamic  roll  angles  of  22.5,  45, 

67.5,  and  90  deg,  in  addition  to  the  0  deg  configuration  shown  in  figure  12.  Several  cases  with 
part  of  the  WT  sting  modeled  were  investigated  to  determine  if  there  were  sting  effects  present. 
The  WT  walls  were  not  modeled.  All  simulations  were  performed  with  the  actuator  wing  at  a 
fixed,  fully  deployed  orientation. 
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Body  0  (no  sting)  Body  1  (no  sting)  Body-Wing-8°  (no  sting) 

^  <=3* 


Body  1  (+  sting)  Body-Wing-8°  (+  sting) 


Figure  12.  Representative  CFD  model  configurations. 

The  computational  domain  was  a  sphere  with  an  outer  boundary  located  50  cal.  from  the 
projectile.  Figure  13  shows  some  views  of  the  mesh  on  the  symmetry  plane  of  the  computational 
domain  and  the  projectile  surfaces.  The  mesh  was  made  using  MIME  version  4.1,  (35)  an 
unstructured  mesh  generator  using  tetrahedral,  prism,  and  pyramidal  cells.  The  mesh  sizes  were 
9.6,  1 1.5,  and  19.7  million  cells  for  the  solid  body  alone,  body  alone  with  empty  actuator  cavity, 
and  body  with  the  actuator,  respectively.  The  results  show  some  effect  of  the  sting  on  the  axial 
force.  Therefore,  the  validation  simulations  were  performed  with  the  sting  (and  at  1.75  scale), 
while  the  concept  characterization  simulations  were  performed  with  the  unsealed  projectile 
model  with  no  sting. 
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Figure  13.  Computational  mesh  on  symmetry  plane  and  projectile  surfaces. 


Validation  simulations  were  performed  at  -12  deg  <  a<  12  deg  in  2-deg  increments  at  Mach  0.2 
at  standard  temperature  and  pressure  (288.15  K  and  101.325  kPa,  respectively).  The  data  from 
the  discrete  angles-of-attack  were  then  least-squares  fit  with  an  appropriate  polynomial 
expansion  in  angle-of-attack.  All  simulations  were  steady-state.  Figure  14  shows  some  typical 
flowfield  visualization  results,  including  the  effect  of  the  sting  on  the  base  flow.  The  wing 
actuator  has  a  pronounced  effect  on  the  base  flow  and  the  sting  has  a  small  effect  on  the  base 
flow.  Figure  15  shows  a  close-up  view  of  the  actuator  region  for  the  four  actuator  wing 
deflections  considered.  Qualitatively,  there  is  not  much  difference  in  the  flowfields,  but  as  shown 
in  figure  10,  the  normal  force  scales  with  wing  deflection.  It  must  be  noted  that  these  flowfields 
are  forced  to  a  steady-state  by  the  solution  technique.  These  flowfields  are  most  likely  highly 
unsteady,  but  since  the  rotation  of  the  actuator  in  this  stage  is  not  included,  the  steady-state 
solution  is  adequate. 
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Figure  14.  Mach  number  contours  on  vertical  symmetry  plane  at  a  =  0  deg. 


Figure  15.  Close-up  view  of  Mach  number  contours  on  vertical 

symmetry  plane  at  a  =  0  deg,  actuator  deflections  of  4,  6,  8, 
and  1 0  deg. 
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Figure  16  shows  the  comparison  of  CFD  with  WT  data  for  the  body-alone  configuration  with  the 
empty  actuator  cavity.  Figure  16a  shows  the  total  axial  force,  for  which  the  data  compare 
reasonably  well  at  a  =  0  deg.  The  case  without  the  sting  shows  a  larger  increase  in  axial  force 
than  observed  in  the  WT  data  as  the  projectile  deflection  increases  or  decreases.  The  case  with 
the  sting,  however,  compares  very  well  with  the  WT  data  at  positive  angles-of-attack  and  much 
better  than  the  case  without  the  sting  at  negative  angles.  Figure  16b  shows  the  forebody  axial 
force,  i.e.,  the  axial  force  on  all  surfaces  except  the  projectile  base.  The  two  CFD  results  nearly 
overlay,  indicating  the  differences  are  all  on  the  projectile  base  and  due  to  the  sting  effects. 
Figures  16b  and  16c  present  the  normal  force  and  pitching  moment  results,  respectively,  and 
indicate  that  there  are  no  sting  effects  on  these  parameters.  The  normal  force  compares  very  well 
to  the  experimental  data  while  the  predicted  pitching  moment  slope  is  slightly  higher  than 
observed  in  the  WT.  The  latter  can  possibly  be  due  to  a  laminar-to-turbulent  transition  in  the  WT 
that  is  not  accounted  for  in  the  CFD,  which  is  assumed  completely  turbulent. 


Figure  16.  Comparison  of  CFD  and  WT  data  for  body  alone  (configuration  1):  (a)  axial  force, 
(b)  forebody  axial  force,  (c)  normal  force,  and  (d)  pitching  moment. 
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Figure  17  shows  the  comparison  of  CFD  with  WT  data  for  the  11120  configuration.  From  the 
axial  force  plots  (figures  17a  and  17b),  it  can  be  seen  that  sting  effects  are  only  observed  at 
a>  2  deg.  Also,  some  sting  effects  can  now  be  observed  on  the  pitching  moment  (figure  16d)  for 
a>  5  deg  as  indicated  by  the  separation  of  the  two  CFD  results.  The  normal  force  (figure  16c)  is 
also  over  predicted  for  a  >  -4  deg.  The  experimental  uncertainty  and  comparison  with  CFD 
indicates  that  the  CFD  results  can  be  applied  with  confidence  for  a  comprehensive 
characterization  of  the  concept. 


Figure  17.  Comparison  of  CFD  and  WT  data  for  configuration  11120:  (a)  axial  force,  (b)  forebody  axial  force, 
(c)  normal  force,  and  (d)  pitching  moment. 

The  purpose  of  the  concept  assessment  phase  of  the  CFD  study  was  to  determine  the  dynamic 
derivatives,  i.e.,  Magnus,  roll-damping,  and  pitch-damping  moments  on  the  down-selected 
configuration  (1 1 120).  The  scale  of  the  projectile  in  this  phase  of  the  study  was  40  mm,  so  the 
mesh  was  scaled  down  appropriately.  Roll-damping  and  Magnus  moments  were  generated  by 
performing  unsteady  simulations  of  the  projectile  rolling  at  the  nominal  spin  rate  (54.2  Hz)  at 
angles-of-attack  -12  deg  <  a<  12  deg.  Pitch  damping  was  calculated  using  both  the  steady-state 
lunar  coning  and  transient  planar  pitching  methods  (36-38). 
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The  Magnus  moment  for  the  body-alone  configuration  was  very  nonlinear,  ranging  from  -0.77  at 
a  =  1  deg  to  +0. 10  at  a  =  12  deg.  The  Magnus  moment  for  the  1 1 120  configuration  was  positive 
and  nearly  constant,  with  values  ranging  from  0.25  to  0.5.  The  pitch-damping  moment  was  more 
difficult  to  accurately  predict  than  other  projectiles  previously  analyzed.  The  combination  of 
nearly  incompressible  flow  and  very  small  pitch-damping  values  (on  the  order  of  1)  due  to  the 
low  length-to-diameter  (L/D)  ratio  is  the  likely  reason.  It  was  difficult  to  obtain  linear  results 
using  the  lunar  coning  method — an  inherent  assumption  of  this  method.  For  the  planar  pitching, 
the  results  were  converged  to  a  lower  tolerance  at  each  time  step  than  past  projectile  geometries 
investigated  at  higher  subsonic  and  supersonic  Mach  numbers.  These  issues  are  discussed  in 
more  depth  in  reference  3 1 .  The  pitch-damping  moments  determined  from  the  transient  planar 
pitching  method  was  more  accurate  than  those  from  the  steady-state  lunar  coning  method.  The 
pitch-damping  values  provided  as  input  to  the  aerodynamic  model  were  -1.24  for  the  body-alone 
configuration  and  -3.04  for  the  11120  configuration. 

All  aerodynamic  techniques  (SEAP,  WT,  CFD)  were  utilized  to  provide  the  aerodynamic 
coefficient  database.  Subject  matter  expertise  was  applied  to  reconcile  the  best  aerodynamic  data 
available  from  various  sources.  The  projectile  body  axial  force  coefficient,  showing  the  WT, 
CFD,  and  merged  (final  database)  data  is  shown  in  figure  18.  Zero-yaw  axial  force  is  similar 
between  WT  and  CFD  methods  and  the  yaw-dependent  axial  force  is  slightly  higher  for  the  CFD. 
Figure  19  provides  the  normal  force  of  the  control.  The  shape  of  the  curves  is  similar  for  both  the 
WT  and  CFD  data.  The  CFD  predicts  about  a  0.05  higher  trim  force  coefficient.  A  summary  of 
the  aerodynamic  database  is  given  in  tables  5-7.  SEAP,  WT,  and  CFD  data  are  compiled  along 
with  the  final  merged  dataset.  Static  aerodynamics,  roll  damping  moment  coefficient,  and  pitch 
damping  moment  coefficient  for  the  projectile  are  in  table  5.  The  nonlinear  Magnus  force  and 
moment  of  the  projectile  are  tabulated  in  table  6.  Table  7  provides  the  aerodynamic  coefficients 
for  the  control  as  a  function  of  the  stowed  roll  angle. 
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Figure  18.  Axial  force  coefficient  for  projectile  body. 


Figure  19.  Normal  force  coefficient  for  control. 


Table  5.  Projectile  aerodynamics. 


Method 

C*o 

Cx  2 

a L 

cNa 

Cn  3 

or5 

Cma 

r 

ma  3 

Ci 

Lp 

r 

^mq 

SEAP 

0.19820 

2.1 

2.08342 

0 

2.09488 

0 

-0.03055 

-3.3 

WT 

0.169782 

1.011851 

1.829601 

-4.94869 

1.932517 

-1.99989 

NA 

NA 

CFD 

0.172793 

3.958777 

2.076615 

-3.07324 

2.065624 

0.423832 

-0.02835 

-1.24 

Merged 

0.17204 

3.22205 

2.01486 

-3.54210 

1.99907 

-0.78803 

-0.02937 

-1.75500 
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Table  6.  Projectile  magnus  aerodynamics. 


Method 

a  (deg) 

0 

1 

2 

3 

4 

5 

7.5 

10 

15 

20 

30 

SEAP 

Cv 

Ypa 

-0.52 

-0.52 

-0.52 

-0.52 

-0.52 

-0.52 

-0.52 

-0.52 

-0.52 

-0.52 

-0.52 

cn 

npa 

-1.6 

-1.54 

-1.37 

-1.13 

-0.88 

-0.7 

0.05 

0.03 

0.03 

0.02 

0.01 

CFD 

Cy 

Ypa 

NA 

0.392990 

225 

0.162194 

013 

NA 

NA 

0.374164 

713 

0.636513 

499a 

NA 

0.860910 

632a 

NA 

NA 

Cn 

npa 

NA 

0.598981 

682 

0.453168 

691 

NA 

NA 

0.124272 

894 

0.014947 
65  8a 

NA 

0.095972 

186b 

NA 

NA 

Merged 

Cy 

rpa 

0.33784 

0.16474 

-0.00835 

-0.14244 

-0.27653 

-0.41062 

-0.57459 

-0.69153 

-0.77568 

-0.77568 

-0.77568 

Cn 

rLpa 

-0.95860 

-0.83424 

-0.68238 

-0.54015 

-0.39543 

-0.26820 

0.00631 

0.04909 

0.07948 

0.07698 

0.07448 

aCFD  data  collected  at  a  =  8  deg. 
bCFD  data  collected  and  a  =  12  deg. 


Table  7.  Control  aerodynamics  for  all  stowed  roll  angles. 


Method 

0, 

CXo 

CX 

Cx  2 

Cx  3 

cr3 

Cn0 

CNa 

Cn  3 

or3 

r 

u m0 

r 

U  ma 

r 

0  deg 

0.10598 

0.48966 

0.60734 

0.08227 

0.20748 

0.46091 

-1.19402 

-0.11250 

-0.20850 

-0.20850 

WT 

30  deg 

0.08936 

0.48966 

0.08423 

0.08227 

0.18315 

0.45047 

1.00463 

-0.12169 

-0.44139 

0.49181 

60  deg 

0.04008 

0.48966 

-0.16194 

0.08227 

0.09635 

0.45974 

4.04593 

-0.08091 

-0.39343 

0.23514 

90  deg 

0.03674 

0.48966 

-0.41301 

0.08227 

0.06707 

0.31557 

8.50209 

-0.06287 

-0.57775 

0.75245 

CFD 

0  deg 

0.11796 

0.65287 

0.78644 

0.10969 

0.21929 

0.46867 

-1.33984 

-0.13927 

-0.30858 

-2.12313 

0  deg 

0.10598 

0.48966 

0.60734 

0.08227 

0.20748 

0.46091 

-1.19402 

-0.11250 

-0.20850 

-0.20850 

Merged 

30  deg 

0.08936 

0.48966 

0.08423 

0.08227 

0.18315 

0.45047 

1.00463 

-0.12169 

-0.44139 

0.49181 

60  deg 

0.04008 

0.48966 

-0.16194 

0.08227 

0.09635 

0.45974 

4.04593 

-0.08091 

-0.39343 

0.23514 

90  deg 

0.03674 

0.48966 

-0.41301 

0.08227 

0.06707 

0.31557 

8.50209 

-0.06287 

-0.57775 

0.75245 

4.  Aerodynamic  Model 


The  aerodynamic  model  relates  the  aerodynamic  coefficient  data  collected  in  experiments  and 
computations  to  the  aerodynamic  forces  and  moments  acting  on  a  projectile  in  flight.  Total 
aerodynamic  forces  and  moments  are  separated  into  rigid  and  moveable  surfaces. 


X  =  XR  +XM 

Y  =  Y R  +  Y m 

(2) 

z  =  ZR  +ZM 

l  =  lr  +  lm 

M  =  Mr  +  Mm 

(3) 

N  =  Nr  +  Nm 

Rigid  aerodynamic  surface  forces  include  static  (linear  and  nonlinear)  and  dynamic  terms  for  the 
projectile  body.  Symbols  in  parenthesis  indicate  functional  dependence  of  aerodynamic 

1  o  71  n 

coefficients.  The  dynamic  pressure  is  q  =  ~pV  and  aerodynamic  reference  area  is  5  =  -D 
where  D  is  the  projectile  diameter  and  V  is  the  total  velocity. 

XR  =  -qs[cXo(M)  +  CXa2(M)a2] 

_  r  o  pD~\ 

Yr  =  -qS  [CNam p  +  CNa3mP3  -  Cv(M,a)a|_]  (4) 

_  r  ,  pD 

ZR  =  -qS  [CNa(M)a  +  (^(Mja3  +  CYpa(M,a W~ 

Rigid  aerodynamic  surface  moments  include  static  (linear  and  nonlinear)  and  dynamic  terms. 

D  T  p£>  1 

L*  =  qSD 

„  r  ^  qD  pD~\ 

MR  =  qS \Cmama  +  Cm^  (M)a3  +  Cmq(M )  ^  +  CUpa(M,  a)  (3  (5) 

„  0  rD  vD~ 

Nr  =  qS  \-Cma{M)(3  -  C^tM)/?3  +  Cmq(M)  —  +  Cnpa(M,a)a *— 

The  following  algorithm  may  be  used  to  calculate  moveable  aerodynamic  surface  forces  and 
moments  for  the  ith  control  mechanism.  First,  relevant  roll  angles  for  this  problem  are  defined. 
Neglecting  greater  than  6  degrees-of- freedom  (DOF)  dynamics,  the  relative  roll  angle  of  the 
control  mechanism  is  <t>i  =  —0  +  4>Cmd  (where  <PCmd  is  the  commanded  roll  direction  of 
maneuver)  and  the  Earth- fixed  roll  angle  of  the  control  mechanism  is  0C  =  0  +  <!>,.  The  stowed 
roll  angle  of  the  ith  control  mechanism  then  becomes  =  <t>j  —  <pB.  (where  0B.  is  the  roll 
angle  of  the  control  mechanism  mounted  in  the  projectile  body). 
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Compute  local  velocity  at  each  moveable  aerodynamic  surface  from  center-of-pressure  (CP)  data 
(xCP  is  axial  CP  and  rCP  is  radial  CP),  geometry  (0B.),  and  6  DOF  states. 

VmJi  =  Vcg/i  +  we/;  x  ^cg^cPj  (6) 


where  VCG/j  =  [u  v  w]T,  a. )B/,  =  [P  q  r]T ?  and 

T 

rcc^cPi  =  1*cp  rCPcos(ct)Bi )  rCPsin((pB.)]  . 


Obtain  velocity  in  ith  moveable  aerodynamic  surface  coordinates  ([uMj  vm-l  wMt]T)  using  the 
transformation  matrix. 


T’eMi  — 


■1  0 
0  cosCdh) 
.0  — sin(<t>  i) 


0 

sin(< J>j) 
cos  (<f>i) 


(7) 


Calculate  total  and  local  moveable  aerodynamic  surface  angle-of-attack. 


&]\/j.  —  CLSITI 


=  asin 


fvMi2  +  wMi2 

UMt 

-V 

2  +  ^Mi2  +  WM;2 

/ 


vMi 


\ 


\luMiz  +  vM/  +  wM.y 


(8) 


Determine  moveable  aerodynamic  surface  axial  force,  normal  force,  roll  moment,  and  pitching 
moment  coefficients.  Side  forces  and  side  moments  on  the  moveable  aerodynamic  surface  are 
neglected. 

,$i)  +  ClM(M,$t)pMl  + 

Cf  =  +  ClK{M,$t)PMl  + 

C  =  C,"(M,fc>  +  («.&)«*, 

c"‘  =  +  CaM(M,  <?,)/?„,  + 

Compute  moveable  aerodynamic  surface  axial  force,  normal  force,  roll  moment,  and  pitching 
moment. 


_|_  QM  f  A/f  £ 


^3m 


AT; 


*Mi  =  -<7MtSC"‘ 
^  =  -qMiSCf 
=  - qM.SDC 
NMi  =  -qMiSDCl 


(10) 
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Transform  moveable  aerodynamic  surface  forces  and  moments  to  body  coordinates. 


\xMi- 

\XMi 1 

YMi 

—  f~l 

~  *  BMt 

YMi 

-ZMi- 

B 

-  0  - 

\XM1] 

YMi 

—  f-1 

—  1BMi 

0 

-ZMi- 

B 

-NMi- 

(11) 


5.  Flight  Mechanics 


The  equations  of  motion  are  formulated  in  the  inertial  frame  while  the  body  frame  was  used  for 
collection  of  aerodynamic  data.  In  this  report,  the  Earth  coordinate  system  (subscript  E)  is  used 
for  the  inertial  frame  and  the  body- fixed  coordinate  system  (subscript  B)  is  used  for  the  body 
frame.  These  right-hand  coordinate  systems  are  related  by  the  Euler  angles  for  roll  (0),  pitch  (0), 
and  yaw  (0)  as  shown  in  figure  20. 
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Rotating  through  the  sequence  of  Euler  angles  yields  the  transformation  matrix  from  quantities  in 
body-fixed  coordinates  to  Earth  coordinates. 


Tbe 


CqC ip  SqSqC^  C^Sfjj  CpSgC.p  +  SfpS^p 

CeSxp  SqSqS^p  +  CfC^p  C^SgC^p  + 

—Sq  S^Cq  CfpCQ 


(12) 


Projectile  flight  states  are  center-of-gravity  position  [x  y  Z]T ,  attitude  [(p  9  ip]T ,  body 

translational  velocity  [u  v  w]T,  and  body  rotational  velocity  [P  <7  r]T .  Translational  and 
rotational  kinematics  for  the  body- fixed  coordinate  system  are  provided  below  (i,  2). 


X 

CqCxI)  S^SqC^  —  C^S^p  C^SgC-^  +  SfpS^ 

-u- 

y 

— 

Gq Sip  S^SgS^  T  Cfp Sq C-^p  + 

V 

-Z- 

—  Sq  SqCq  C^Cg 

-w- 

0 

1  S^tg 

c<j>tg 

-p- 

6 

= 

0  C0 

q 

A 

0  S<p/Cg 

Cep/  Cg 

-r- 

(13) 


(14) 


Newtonian  kinetics  may  be  applied  to  a  projectile  in  free-flight.  The  translational  and  rotational 
dynamics  are  expressed  in  body-fixed  coordinates  ( 1 ,  2). 


u 

v 

w. 


\X  +  XGl 

-  0 

—r 

q ' 

-u- 

Y  +  Yc 

— 

r 

0 

-p 

V 

z  +  ZG_ 

-q 

P 

o . 

-w- 

V 

L' 

1 

o 

-p- 

q 

=  r1 

M 

-I-1 

r  0  — p 

7 

q 

-f. 

-N. 

—q  p  0  . 

-r- 

(15) 


(16) 


The  forcing  functions  in  the  governing  equations  have  components  due  to  aerodynamics  and 
gravity.  The  aerodynamic  model  is  provided  earlier.  Gravity  acts  only  as  a  force 
([XG  Yg  ZgY  =  Ttbe[  0  0  gf). 


6.  Delivery  in  the  Ballistic  Environment 


The  process  of  delivering  a  lethal  payload  via  a  man  portable  weapon  in  the  ballistic 
environment  must  be  described  and  modeled.  First,  the  target  location  (and  sometimes  gun 
location)  is  measured  or  estimated.  For  direct  fire  engagements  (line  of  sight),  these  data  are 
obtained  from  a  rangefinder  device  or  a  human  guess.  Targeting  information  comes  through  a 
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network  from  a  collective  asset  for  indirect  fire  (or  more  complex  direct  fire)  engagements. 
Regardless  of  the  targeting  method,  there  is  uncertainty  associated  with  the  data. 

A  fire  control  solution,  or  desired  weapon  pointing  angle,  is  determined  based  on  the  targeting 
data.  The  weapon  is  aimed,  sometimes  through  the  use  of  an  iron  sight  or  an  aim-assist  device, 
and  the  projectile  launched.  Errors  are  associated  with  these  processes  because  the  fire  control 
solution  is  never  perfect,  the  weapon  cannot  precisely  achieve  the  desired  aiming  angles,  and  the 
launch  has  variability  in  the  muzzle  velocity  and  angular  rates  of  the  projectile. 

The  projectile  flies  through  the  atmosphere  on  the  way  to  the  target.  Atmospheric  density,  sound 
speed,  and  wind  vary  throughout  time  and  space  and  influence  projectile  flight.  Another  factor 
that  affects  the  projectile  flight  is  physical  variability.  The  mass  properties  and  aerodynamics 
change  from  round  to  round  due  to  effects  such  as  manufacturing  tolerances. 

The  fuze  prescribes  the  terminating  location  of  the  projectile.  A  variety  of  fuzing  mechanisms, 
point-detonating,  proximity,  time,  or  tums-count,  are  available.  The  method  and  uncertainty 
dictate  the  burst  point  of  the  warhead. 

Table  8  provides  the  models  and  input  data  used  in  this  paper  for  these  delivery  processes. 
Uniform  (U)  and  normal  (M)  distributions  are  used  throughout  the  modeling.  The  targeting 
module  assumes  a  rangefinder  device  estimates  the  range-to-target  ( RT )  which  in  turn  corrupts 
the  fire  control  solution.  Aiming  assumes  a  supported  weapon.  Spin-stabilized  projectile  drift  is 
accounted  for  in  the  azimuth  ( 1 ,  2).  Muzzle  velocity  is  modeled  based  on  recoil  impulse 
limitations  and  propulsion  efficiency.  The  spin  rate  at  launch  is  prescribed  by  the  muzzle 
velocity  and  gun  twist  to  within  variation  of  the  rifling  engagement.  Atmospheric  density  is 
perturbed  from  the  1962  International  Standard  Atmosphere  (ISA).  Wind  is  modeled  as  a 
constant  horizontal  component  over  the  flight  in  any  direction.  All  fuzes  burst  at  a  nominal  point 
for  that  fire  control  solution  (e.g.,  range-to-go  [/?50]  for  proximity,  time  to  range  [tf]  for  time, 
roll  angle  to  range  [<pf]  for  turns)  to  within  some  tolerance  due  to  relevant  effects  such  as  radio¬ 
frequency  ranging,  clock  errors,  or  turn  quantization. 
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Table  8.  Models  and  data  for  man  portable  delivery  process. 


Category 

Process 

Model 

Data 

Targeting 

Targeting 

Rt  =  U (Rt>  i Gr) 

ctr  —3m 

Launch 

Weapon  elevation 

90  =  N(fiRTla0) 

table  90  =  f(RT ),  oe  =  0.0014  rad. 

Weapon  azimuth 

V'o  = 

table  xpQ  =  f(RT,  90 ),  a \p  =  0.0014  rad. 

Velocity 

<0 

41 

ii 

IR  =  24.5  N-s,  fP  =  1.1,  gv  =  2  m/s 

Spin  rate 

2n  2n 

H  =  30  calibers/rev,  op  =  0.0035  % 

Pitch  rate 

o 

II 

O 

£ 

<Jq=  2  rad/s 

Yaw  rate 

r0  =  FJ(0,  <rr) 

or  =  2  rad/s 

Atmosphere 

Density 

p(z)  =  N(1962  ISA,  Op) 

op  %  0.015% 

Wind  speed 

Vw=  N(0  ,aVw) 

ov  =3.36  m/s 

vw 

Wind  direction 

ew  =  u(o,ct0w) 

Of)  =  2n  rad 

Mass  Properties 

Mass 

m  =  N(m,  om) 

in  =  0.342  kg,  om  =  0.5  % 

Diameter 

D  =  N(D,  aD) 

D  =  0.040  m,  oD  =  0.1  % 

Axial  inertia 

Ix  =  N(Ix,<7jx) 

Vx  =  8.54e-5  kg- nr,  oIx  =  3  % 

Transverse  inertia 

E~h 

<4; 

& 

II 

Ft  =  3.56e-4  kg-m2,  o,T  =  1  % 

Aerodynamics 

Axial  force 

Cx  —  N(  Cx,aCx) 

Cx  (see  aero  tables),  oCx  =  2  % 

Normal  force 

S: 

la 

(<J 

II 

<o 

CN  (see  aero  tables),  oCn  =  5  % 

Pitching  moment 

Cm  =  N(Cm,C7Cm) 

Cm  (see  aero  tables),  oCm  =  4  % 

Roll  damping  Moment 

Clp  =  N(CTp,aQp) 

Cip  (see  aero  tables),  oc  =  5  % 

Pitch  damping 
moment 

r 

-0 

II 

p) 

> 

Cmq  (see  aero  tables),  crCm^  =  15  % 

Magnus  force 

r 

R 

II 

/ — s 

«? 

R 

C^a  (see  aero  tables),  =  ^  % 

Magnus  moment 

Z*— S 

a 

y 

II 

a 

a 

vf 

C^a  (see  aero  tables),  o cnva  =  15  % 

Fuzing 

Proximity 

Rgo  ^(Rgo’^Rgo) 

Rgo  =  3  m,  aRg0  =  0.5  m 

Time 

tf= 

table  if  =  f{RT ),  ot  =  0.05  % 

Turns-count 

4>f  =  U(/(Rr),±a0/) 

table  (j)f  =  f(RT ),  o^f  =n/2  rad 

7.  Delivery  Characterization 


The  ballistic  delivery  accuracy  was  assessed.  The  models  and  input  data  previously  described  for 
the  targeting,  launch,  atmosphere,  mass  properties,  aerodynamics,  and  fuzing  were  implemented 
in  a  Monte  Carlo  simulation  of  the  flight  mechanics.  The  target  was  placed  at  200  m.  Burst 
points  were  compiled  with  point-detonating,  proximity,  time,  and  tums-count  fuzes.  Figures  21-22 
show  three-dimensional  burst  points  for  ballistic  flight.  All  point-detonating  fuze  cases  impact  on 
the  ground  plane.  Time  and  tums-count  fuzes  stand  the  bursts  up  more  in  the  vertical  plane.  All 
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results  feature  significant  spread  in  the  downrange  direction  (usually  more  than  10  m)  since  the 
launch  is  relatively  low  elevation  (approximately  15  deg).  The  crossrange  dispersion  was  usually 
less  than  1-2  m. 
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Figure  21.  Horizontal  plane  ballistic  burst  points. 
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Figure  22.  Vertical  plane  ballistic  burst  points. 
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Ballistic  delivery  accuracy  provides  a  metric  for  initial  assessment  of  the  feasibility  of  this  small 
diameter  guided  projectile  concept.  The  control  authority  afforded  by  this  concept  must  be  larger 
than  the  ballistic  dispersion.  The  aerodynamic  data  collected  in  the  WT  (contained  in  tables  3-4) 
was  used  to  quantify  the  maneuver  footprint  of  the  various  configurations.  Flight  simulations 
were  conducted  with  nominal  launch,  atmosphere,  mass  properties,  and  aerodynamics  with  an 
elevation  that  reaches  a  ballistic  range  of  200  m.  Maneuvers  began  at  launch  and  continued 
throughout  flight  to  ground  impact.  Maneuver  direction  was  varied  across  eight,  equally 
distributed  angles  to  provide  a  sense  of  crossrange  and  downrange  course  corrections  available. 


Figure  23  demonstrates  the  control  authority  footprint  for  the  different  airfoil  profiles  with  a  two- 
actuator  variant.  In  general,  the  projectile  flies  furthest  downrange  for  a  maneuver  up  since  this 
represents  the  longest  time-of-flight.  The  difference  between  the  maneuver  down  and  the 
ballistic  impact  point  is  larger  than  the  differences  for  the  maneuver  left  or  right  mainly  due  to 
the  angle-of-fall  (near  20  deg  from  horizontal).  The  greater  normal  force  provided  by  the 
cambered  airfoil  yields  significantly  more  control  authority.  Comparing  the  horizontal  plane 
control  authority  for  the  NACA  4412  airfoil  with  the  ballistic  burst  points  suggests  that 
maneuverability  is  sufficient  to  account  for  sources  of  delivery  error. 
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Figure  23.  Ground  plane  control  authority  for  airfoil  parameter. 
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The  total  downrange  and  crossrange  course  corrections  were  calculated  from  the  horizontal  plane 
maneuver  impacts  for  the  airfoil  parameter.  These  data,  provided  in  figure  24,  quantify  an  almost 
factor  of  2  increase  in  the  downrange  and  crossrange  control  authority  for  the  cambered  airfoil 
over  the  symmetric  airfoil. 
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Figure  24.  Control  authority  for  airfoil  parameter. 

The  downrange  and  crossrange  maneuver  authority  was  calculated  for  a  wider  parameter  set  in  a 
manner  similar  to  that  shown  for  the  airfoil  parameter.  Figure  25  illustrates  that  the  larger 
planform  of  the  wing  increases  the  control  authority.  Examination  of  the  results  in  figure  26 
indicates  that  control  authority  increases  nearly  linearly  with  the  wing  deflection  angle.  The  data 
for  the  standoff  distance  shown  in  figure  27  quantifies  how  getting  the  wing  further  away  from 
the  projectile  increases  maneuverability.  Finally,  adding  more  wings  as  provided  in  figure  28 
linearly  increases  the  control  authority. 
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Figure  25.  Control  authority  for  planform  parameter. 
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Figure  26.  Control  authority  for  deflection  parameter. 
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Figure  27.  Control  authority  for  standoff  parameter. 
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Control  authority  was  larger  than  the  ballistic  dispersion  for  a  significant  subset  of  the  wing 
parameters.  The  results  of  the  parametric  study  of  maneuverability  and  practical  packaging 
considerations  suggest  a  best  set  of  wing  characteristics.  These  parameters  included  a  two- 
actuator  variant  with  a  NACA  4412  airfoil  sized  to  0.48-cal.  chord  x  0.61 -cal.  span  deflected  to 
8  deg  with  maximum  standoff.  The  CFD  aerodynamic  predictions  were  performed  with  this 
configuration  and  the  results  presented  in  the  remainder  of  this  report  possess  the  aerodynamic 
coefficients  for  this  geometry  (as  given  in  tables  5-7). 

Typical  maneuver  flight  response  is  outlined  in  figures  29-33.  The  down-selected  configuration 
was  launched  at  approximately  15  deg  under  nominal  conditions  and  an  up  maneuver  was 
commanded  throughout  the  flight.  The  pitch  angle-of-attack  history  shows  that  the  maneuver 
excites  the  precession  and  nutation  modes.  Precession  is  evident  in  the  lower  frequency,  higher 
amplitude  oscillation  angular  motion.  The  maneuver  concept  induces  an  angle-of-attack  in  the 
projectile  just  over  4  deg.  The  yaw  angle-of-attack  demonstrates  similar  precession  and  nutation 
motion.  Yaw  of  repose  over  about  1  deg  for  this  spin  rate  and  launch  elevation  can  be  seen.  This 
projectile  features  a  launch  spin  rate  of  about  54  Hz  (corresponding  to  a  gyroscopic  stability 
factor  around  2.3)  and  decreases  about  2  Hz  over  4  s.  Mach  number  decreases  throughout  flight 
from  just  under  0.2  at  launch.  The  vertical  plane  trajectory  demonstrates  that  the  range  is 
extended  over  50  m  past  the  ballistic  impact  point  (200-m  ballistic  range  at  this  elevation)  from 
the  maneuver  concept  inducing  a  4-deg  angle-of-attack.  No  deleterious  flight  instabilities  were 
encountered  during  flight  simulations. 


Figure  29.  Pitch  angle-of-attack  for  maneuver  up. 
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Figure  30.  Yaw  angle-of-attack  for  maneuver  up. 


Figure  3 1 .  Spin  rate  for  maneuver  up. 
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Figure  33.  Vertical  plane  trajectory  for  maneuver  up. 
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8.  Conclusions 


This  study  demonstrated  the  utility  of  a  new  concept  for  maneuvering  a  small-diameter 
projectile.  The  concept  was  detailed  and  the  aerodynamics  of  this  shape  was  determined  using  a 
variety  of  techniques.  WT  experiments  were  conducted  on  numerous  wing  geometries  for  initial 
feasibility.  Experimental  data  were  used  to  down-select  the  configuration  and  validate  CFD 
predictions  prior  to  applying  the  computations  for  more  advanced  purposes  such  as  determining 
dynamic  derivatives. 

An  aerodynamic  model  to  relate  WT-  and  CFD-derived  coefficients  to  free-flight  forces  and 
moments  was  derived  specifically  for  this  novel  configuration.  Effects  such  as  the  time- 
dependent  (or  roll  history)  exposure  of  the  wing  to  the  airstream  were  included  in  the  model. 

The  nonlinear  mechanics  underpinning  flight  simulations  were  presented.  The  delivery  process 
and  modeling  of  a  lethal  payload  via  this  guided  flight  concept  was  summarized.  Physics-based 
models  and  input  data  for  targeting,  and  variation  due  to  launch,  atmosphere,  physical 
characteristics,  and  fuzing  were  outlined.  Delivery  characterizations  were  performed  via 
simulation  to  assess  control  authority  requirements. 

Realistic  maneuver  footprints  were  obtained  through  flight  simulation.  Trade  studies  were 
conducted  to  quantify  the  relationship  between  wing  parameters  and  control  authority.  A 
characteristic  set  of  flight  dynamic  metrics  were  provided  for  understanding  guided  flight 
behavior. 

This  study  indicates  that  control  authority  of  this  novel  concept  is  sufficient  to  compensate  for 
ballistic  delivery  error  sources.  These  results  were  obtained  despite  the  challenges  associated 
with  flying  a  small-diameter,  spin-stabilized  projectile  with  low-dynamic  pressure.  The 
components  involved  with  implementing  this  concept  are  low  cost  (e.g.,  commercial  DC 
motors).  These  technologies  support  realization  of  a  new  class  of  precision  weapons  for  the 
Army  at  the  squad  level. 

Future  computational  and  experimental  efforts  focus  on  understanding  the  effects  of  the  transient 
wing  exposure  to  the  airstream  on  the  flight  behavior.  Additionally,  mechatronic  design  and 
evaluation  (to  include  gun  hardening)  of  the  actuation  technology  must  be  performed. 
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AFRL 

ARL 

CFD 

CP 

DOD 

DOF 

DSRC 

ERDC 

HLLC 

ISA 

L/D 

NACA 

PRODAS 

RANS 

SEAP 

SST 

TVD 

WT 


U.S.  Air  Force  Research  Laboratory 

U.S.  Army  Research  Laboratory 

computational  fluid  dynamics 

center-of-pressure 

Department  of  Defense 

degree(s)  of  freedom 

DOD  Supercomputing  Resource  Center 

U.S.  Army  Engineering  Research  and  Development  Center 

Harten-Lax-Van  Leer-Contact 

International  Standard  Atmosphere 

length-to-diameter 

National  Advisory  Committee  for  Aeronautics 

Projectile,  Rockets,  and  Ordnance  Design  and  Analysis  System 

Reynolds- Averaged  Navier-Stokes 

semiempirical  aeroprediction 

Shear  Stress  Transport 

T  otal- V  ariation-Diminishing 

wind  tunnel 
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